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Abstract. We consider the Ginzburg-Landau equation, <9 ( u = c^w+w — w|w| , with complex amplitude u(x, t). We 
first analyze the phenomenon of phase slips as a consequence of the local shape of u. We next prove a global theorem 
about evolution from an Eckhaus unstable state, all the way to the limiting stable finite state, for periodic perturbations of 
Eckhaus unstable periodic initial data. Equipped with these results, we proceed to prove the corresponding phenomena 
for the fourth order Swift-Hohenberg equation, of which the Ginzburg-Landau equation is the amplitude approximation. 
This sheds light on how one should deal with local and global aspects of phase slips for this and many other similar 
systems. 
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1. Introduction 

In the last few years we have seen continuous progress in the rigorous understanding of pattern- 
forming systems, describing successfully fronts, modulation equations, nucleation mechanisms 
and scaling phenomena. The aim of the present paper is to add to this body of results a rigorous 
study of "phase slips," that is, of solutions whose number of zeros varies as a function of time. 
The problem of the variation of the number of zeros in parabolic PDE' s has been studied earlier, 
see, e.g. [A], and the references therein. Our paper will deal with two novel aspects of this 
problem, one local and the other global, which use complex order parameters and apply also to 
higher order differential operators. 

The local theory is based on the prototype of the Ginzburg-Landau (GL) equation: 

d t u(x,t) = d 2 u(x, t) + u(x, t) — u(x, t) \u(x, t)\ 2 , (1-1) 

where u : R x R + — >■ C. It allows one to predict from the local shape of the solution at time 
t = whether there will be a phase- slip — a of \u(-, i)\ — at some later time t. The condition 
is that 

u(x,0) \ax 2 + ibx — ad , (1.2) 

where a > 0, and Re b ^ 0. We can then show the appearance of a phase slip for sufficiently 
small d = d(a, b) > 0. (This condition says that the initial condition is "just before" a phase 
slip.) The point that we wish to emphasize is that while we have stated this result as a theorem 
about the initial value problem, in fact it implies that whenever any solution of the GL equation 
gets close to zero it will undergo a phase slip by the mechanism we describe, provided that 
it has a non-zero imaginary part, and that its curvature is small and has a prescribed sign. 
We conjecture that except for very special (i.e., "non-generic," in the appropriate topology) 
solutions, all phase slips occur via the mechanism we describe. To see what we mean by 
"special," note that real valued solutions, for instance, do not satisfy our criteria, but from the 
point of view of considering the GL equation as an amplitude equation, which is the source of 
most of our interest, real solutions are very atypical. The proof of the above result is given 
in Theorem 2.7. It is a perturbative argument, and therefore it can be extended to other, more 
general problems, of which the GL equation is an "amplitude equation." We illustrate this in 
Section 4 for the specific example of the Swift- Hohenberg (SH) equation 

d t U = (Se 2 - (1 + d 2 x ) 2 )U - U 3 . (1.3) 

We consider initial data of the form 

U(x,0) = 2eRe (u(^-ex,0)e ix ^) , (1.4) 

where u satisfies the Ginzburg-Landau slip condition Eq.(1.2). Then, if e > is sufficiently 
small, the solution U(-,t) will have a phase slip, i.e., a pair of roots of U will appear or disappear 
at some time t d/(Se 2 ). This is the way in which the phase slip of the amplitude part of an 
equation translates to the equation itself, see Theorem 4.1. 



Eckhaus Instability 



3 



In many experiments phase slips are produced not by putting the system near the state 
described by Eq.(1.2) but rather by starting with the system in an unstable state, perturbing it 
slightly and waiting for it to undergo a phase slip. This is a global problem, and, obviously 
more care and control are needed for a rigorous treatment than in the local problem. We 
study this global aspect in the setting of the celebrated Eckhaus instability for the GL equation 
in Theorem 3.2. In order to have good control, we consider an Eckhaus unstable stationary 
solution u(x,0) = y/l — q\e q ° x which we perturb by adding a single frequency ee iqiX , or, 
more generally, an initial condition of the form 

^l-q 2 e iq ° x + ^ e n e i(qo+n(qi - qo))x , (1.5) 

with \e n \ small. (Note that the subspace of such functions is left invariant by the GL evolution 
Eq.(l.l).) If q and q x are suitably chosen as described in Section 3, we can show that the 
time-evolution must bring the solution from the initial state to the state y/l — qfe iqiX (up to a 
phase). Thus, this demonstrates the full evolution of an Eckhaus unstable state. In the course of 
this evolution, the wavelength will change from 2n/q to 2n/q 1 and phase slips will occur. 

In Section 4, we sketch, in quite some detail, the proof of the following novel global result: 
Consider the SH equation with initial data 

U(x,0) = 2eRe ^A(^-ex,0)e ix ^j . (1.6) 

Assume the initial data Eq.(l .6) have an amplitude A which is Eckhaus unstable and is perturbed 
as in Eq.(1.5). Then the SH equation will globally change its amplitude modulation as time goes 
on, and come close to a new periodic state whose amplitude has wavevector q 1 . These results 
are based on recent estimates of the relation between the evolution of the SH equation and the 
GL equation [CE2, Sch]. 

2. The Ginzburg-Landau Equation 

In this section we consider the Ginzburg-Landau equation, 

d t u(x,t) = d 2 u(x,t)-\-u(x,t) — u(x,t)\u(x,t)\ 2 , (2-1) 

where u : R x R + — >■ C. Our understanding of the phase-slip mechanism is based on the 
following simple observation. Assume that u(x,t) is a solution satisfying u(Q, t ) = for some 
t > 0. Then, using the well-known regularity properties of solutions of parabolic equations, 
we can expand u near x = 0, t = t as 

u(x,t + tf ) = ^ax 2 + fix + 7T + Odrl 3 / 2 + |a;| 3 ) , 

where a,/3,7 £ C and the term of the form tx has been absorbed in the remainder. Since u 
satisfies Eq.(2.1) and vanishes at x = 0,t = t , we deduce that a = 7. Also, by the phase 
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invariance of the equation, we may assume without loss of generality that a £ R, a > 0. 
Therefore, we expect the function 

u s (x,t) = ^ax 2 + fix + a(i — t ) , a £ R + , /3 £ C , 

to be a good approximation to tf), near a: = 0, < = <„. Note that u s (a;, ^) is an exact solution 
of the linear equation d t u(x,t) = dlu(x,t). Furthermore, if a > and Im/3 ^ (which is 
generically expected to hold), then u 3 (x,t) = if and only if x = and t = t . A sequence of 
plots indicating the phase slip for u s is given in Fig. 1. 



Re u 




Fig. 1: The real and imaginary parts of the function u s (x, t) for 3 successive times. 

We now change our perspective and formulate the problem in terms of initial data. Let 

u (x) = \ax 2 + ibx — ad , (2-2) 

with a > 0, d > 0, Re 6^0. We already know that the corresponding solution 

u Q (x,t) = | ax 2 + ibx + a(t — d) (2-3) 

of the heat equation undergoes a phase slip at x = 0, t = d. The aim of this section is to show 
that this property remains true for any solution u(x,t) of the GL equation with initial condition 

u(x,0) = u (x) + G(\x\ 3 ) , (2.4) 
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near x = 0, provided d is sufficiently small. Our main result can be outlined as follows: 

Proposition 2.1. If d is small, and u(x, 0) is given by Eq.(2.4) near x = 0, then the solution 
u(x,t) oftheGinzburg-LandauequationEq.(2.1)willhaveaphaseslipattimet* = d+(D(d 3 / 2 ) 
and at x* = 0(d 3/2 ). 

Remark. A precise formulation will be given in Theorem 2.7 below. 
Proof. We shall give the proof for the slightly more general equation 

d t u(x,t) = d 2 x u{x,t) + Af(u(x,t)) , (2.5) 

where Af(z) = zP(\z\ 2 ) and P is a real polynomial whose leading coefficient is strictly negative. 
Under these assumptions, every bounded solution u(x,t) of Eq.(2.5) satisfies for sufficiently 
large t: 

\\u(.,t)\\ c3 < U* , (2.6) 

see [CE2]. In the sequel, we always assume that the bound (2.6) holds for all t > 0. Then there 
is a constant N 1 such that 

\Af(z)\ < JVjz| , \d £ z Af(z)\ < for \z\ < U* and 1 = 1,2,3 . (2.7) 

As announced in Eq.(2.4), we consider initial data of the form u(x, 0) = ^(a:) + w (x), 
with ^(,(3;) of the form Eq.(2.2) (on the whole line) and w (x) = C(|a;| 3 ) near x = 0. More 
precisely, we assume that there is a constant W* > such that 

|w (a;)| < W*\x\ 3 , \w' (x)\ < W*\x\ 2 , |<(a;)| < , x e R . (2.8) 

Setting £ = d 2 x and applying the Duhamel formula, we obtain the following representation for 
the solution: 

= (e tC (u 0+ w ))(x) + J*ds (e^ C Af(u(;s)))(x) 

= u (x,t) + iw(a:,tf) + z(a;,^) . 

By construction, u (a;, ^) vanishes at < = <Z and a; = (and only there). Using the Implicit 
Function Theorem, we shall show that the full solution u(x,t) has a unique zero inside the 
region 

Tl d = {(x,t) I x 2 < d, < t < 2d} , (2.10) 

if d is sufficiently small. We bound the three terms in Eq.(2.9) as follows: 

Bounds on u . Since u (x,t) is explicitly given by Eq.(2.3), it is straightforward to prove the 
following 

Lemma 2.2. There is a constant K = K (a, b) such that for d £ [0,1] one has 

K0M)| < K.d 1 ' 2 , |6>^ o 0M)l < K , \d t u (x,t)\ < K , (2.11) 
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forall(x,t) e K d . 

Bounds on w. Here, we exploit the fact that w vanishes at the origin. For any function /, we 
have the explicit representation 



/oo 
dye 
-oo 

If |/(a;)| is bounded by \x\ m , m > 0, we find 

/oo 
dye -(*-v) 2 /M 
-OO 



f(y) 



(2.12) 



y\ 



< c 



< c 



dZe- (2/ *\x + £t 



+ t 1 ' 2 



1/2 I 



(2.13) 



Here and in the sequel, C denotes a constant which can vary from equation to equation. In view 
of (2.8), we can apply Eq.(2.13) to w and its derivatives and obtain 



\w(x,t)\ < C 



\x\ + t 



1/2 



\x\ + t 



1/2 



, \d 2 x w(x,t)\ < C 



\x\ + t 



1/2 



(2.14) 



, \d x w(x,t)\ < C 

Therefore we have shown: 

Lemma 2.3. There is a constant K 1 = K 1 (W*) such that ford £ [0, 1] one has 

\w(x,t)\ < K^ 12 , \d x w(x,t)\ < K x d, \d t w{x,t)\ < K x d 112 , 
forall(x,t) g Tl d . 

Bounds on z. To bound the nonlinear term, we use the integral equation 

z(x,t) = [ ds(e it - s)C Af(u (-,s) + w(-,s) + z(-,s)))(x) , (2.15) 
Jo 

which follows immediately from Eq.(2.9). Note that the right-hand side depends on the solution 
u(x,t) on the whole real line, whereas we are only interested in bounding z for (x,t) £ lZ d . 
We use the following "localization" lemma: 

Lemma 2.4. There is a constant K 2 such that for all f £ C 1 and all t > one has 

\(e iC f)(x)-f(x)\ < K 2 Vt\\f\\ c ,. 



Remark. To be specific, we define 

| cl = max(sup |/(jc)| , sup|/'(a;)|) 
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Proof. Using the representation (2.12) and the Mean Value Theorem, we obtain 



/(*) + 



where i?(a:, £) is some point between a; and a; + £. The assertion follows. 

We apply this lemma to bound the right-hand side of Eq.(2.15). Using Eqs.(2.6), (2.7), 
we find 



e {t - s)c M{u{-, s))j (x) - Af(u(x, s)) < K 2 Vt^l\\Af(u(; s))\\ cl 
< K 2 Vt^N 1 \\u(-,s)\\ cl < K 2 Vt^N 1 U*. 

Integrating this bound over s, we obtain 



/ ds(e it - s)C M{u{-,s))){x)- I dsM{u{x,s)) 
Jo Jo 



< K 3 t 3/2 



(2.16) 



where K 3 = {2/2,)K 2 N 1 U\ 

On the other hand, using Lemma 2.2, Lemma 2.3 and Eqs.(2.7), (2.9), we have for all 
(x,s) G Tl d 

\Af(u(x,s))\ < N^luoix^s)] + \w(x,s)\ + \z(x,s)\) 
< N 1 (K d 1/2 +K 1 d 3/2 + \z(x,s) 



and hence 

/ ds\Af(u(x,s))\ < N^Kq + K 1 )d 1/2 t + N x I ds\z{x,s)\ . (2.17) 
Jo Jo 

Finally, combining Eqs.(2.15)-(2.17), we see that 

\z(x,t)\ < K 4 d 1/2 t + N 1 t ds\z(x,s)\ , 

Jo 

forall(a:,tf) G H d , where K 4 = \[2K 3 +N 1 (K +K 1 ). Therefore, it follows from Gronwall's 
inequality that 

\z(x,t)\ < if 4 tZ 1/2 ^exp(iV 1 ^ < 2K i d 3/2 exp(2N 1 d) , 



for all (x ,t) e n d . 
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To bound the derivatives of z, we differentiate Eq.(2.15) with respect to x and obtain inte- 
gral equations for d x z and d 2 x z. Repeating exactly the same estimates, we find that \d x z(x,t) | < 
Cd and \d?z(x,t)\ < Cd 1/2 for all (a;,*) e K d . Since d t z(x,t) = d 2 x z{x,t) + Af(u(x,t)), we 
see that \d t z(x,t)\ < Cd 1 ^ 2 when (x,t) £ 1Z d . Thus, we have shown: 

Lemma 2.5. There is a constant K 5 = K 5 (a, 6, U* , ) such that for d £ [0,1] one has 

|zOM)| < if 5 ^ 3/2 , |^om)I < # 5< z, |c^om)I < i^ 2 , (2.i8) 

forall(x,t) £ 7£ d . 

We summarize Lemma 2.3 and Lemma 2.5 in the following 

Proposition 2.6. Assume thai: ^(£,0) = u (x) + w (x), with u given by Eq.(2.2) and w 
satisfying Eq.(2. 8). Assume that the solution u(x,t) ofEq.(2.1) satisfies \\u(-,t)\\ c3 < U* for 
all t > 0. Then there is a constant K* = K*(a,b,U* ,W*) so that, for all d £ [0,1], the 
inequalities 

\u(x,t) - (\ax 2 + ibx + a(t - d))\ < K*d 3/2 , 
\d x u(x,t) - {ax +ib)\ < K*d, 

\d t u(x,t) - a\ < K*d 1/2 , 

hold on 1Z d . 

Existence of the phase slip. In view of Proposition 2.6, the existence of a phase slip is now a 
straightforward application of the Implicit Function Theorem. Indeed, we can view the function 
u(x,t) = (Re u(x,t) ,Im u(x,t)) as a (smooth) map from 1Z d to R 2 . By Proposition 2.6, we 
know that ^(0, d) = 0(d 3 ^ 2 ) and the derivative Du(x,t) satisfies 



^ / .\ i ax — Imfe a\ ,i/o> 
Du(x,t)= { Reb Q j + 0{SI 2 ) 



Therefore, since a > and Re b ^ 0,Du(x,t) is close to a non-singular matrix if d is sufficiently 
small. This implies that u is one-to-one from lZ d onto a neighborhood of size O(d) around 
u(0,d) in R 2 . In particular, there is a unique point (x* ,t*) £ 1Z d such that u(x*,t*) = 0. 
Applying the implicit function theorem, we find 

(;>C)-(-<m»)-g r ::^j) + ^». 

Thus, our final result is: 

Theorem 2.7. Assume that u(x,0) = u (x) + w (x), with u given by Eq.(2.2) and w 
satisfying Eq.(2.8). Assume that the solution u(x,t) of the Ginzburg-Landau equation with 
initial data u(x,0) satisfies \\u(-,t)\\ c3 < U* for all t > 0. Then there is a constant d* = 
d*(a,b,U* ,W*) > such that, for all d £ [0,cT], u(x,t) has exactly one phase slip in the 
region K d , . This phase slip will take place at time t* = d+ (D(d 3/2 ) and at x* = + (D(d 3/2 ). 
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3. The Eckhaus Instability 

The considerations of the preceding section were local in nature. Namely, we identified a local 
shape of the amplitude that looks "like just before a phase slip" and then showed that indeed 
a phase slip is going to occur in a well controlled region of space-time. Now, we wish to 
address the more difficult question of global aspects. We will give a mathematical proof 'of a 
phenomenon which is universally believed to hold in a very general context. However, our own 
results here apply in a more restricted setting, due to the requirement of rigor. 

We consider again the Ginzburg-Landau equation, and more specifically perturbations of 
the space periodic time-independent solutions 

s qo (x) = e*°*y/l^q$, ?„G[-1,1]. 

As explained in the Introduction, we shall perturb s qo (x) by adding essentially a single frequency 

ee iqiX , and we shall show that, under suitable assumptions on the parameters q and q 1 , this 
perturbation will grow and eventually drive the system to the new stationary solution s (x). 
Since the wave-lengths 2ir/q , 2w/q 1 are different, phase slips will necessarily occur during 
this transition. 

We first recall that the spectrum of the linearized equation around s (x), considered on 
the subspace of perturbations which are of the form 

u(x) = e iqoX [Ae i6x + 5e" ifl j , A,B e C , 

is given by 

A ± ( 9o ,*) = -(l-? 2 )-* 2 ±\/( 1 -9o 2 ) 2 +4?o 2 * 2 - I 3 - 1 ) 

The quantity A + can take positive values if and only if q% > 1/3, which is the celebrated 
instability region of Eckhaus [CE1]. Therefore, the perturbation ee iqiX will grow at onset if 
we choose q ,q 1 so that q 1 = q — 8 and \ + (q ,8) > 0. We thus take q > l/y/3 and 
< 8 < y/6ql - 2. 

Of course, this linear condition does not guarantee that the system will converge to the 
stationary solution s (x) as t — >■ oo. One reason for this is that Fourier modes other than 
e iqoX , e iqiX will be excited by the nonlinearity, some of them possibly in the Eckhaus stable 
region, and it is not clear a priori which one will eventually "win." To illustrate this point, let 
us fix q = y/2/5, 8 = q /2, and consider the solution u(x,t) of the GL equation with initial 
condition u( x, 0) = s (x) + ee iqiX for some e <C 1. As is easily seen, u(x, t) is of the form 

u( x ^) = ^2 u n( t ) eiqnX > <ln = % -nS , (3.2) 

tiGZ 

where u (0) = y/l — q%, m 1 (0) = e and u n (0) = otherwise. Our choice of q , 8 implies that 
e iqiX is the most unstable Fourier mode around s ( x) (i.e., \ + (q , 8) is maximal for 8 = q /2), 
whereas e iq2X is linearly not excited at all. This is visible in Fig. 2. 
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Fig. 2: The eigenvalue A + (g , 8), describing the instability of a perturbation with wave vector q Q + 8 of the periodic 
stationary solution s . The value q Q = is chosen such that A + is maximal for 8 = ±q /2. It is for 8 = ±q Q . 

Nevertheless, a numerical calculation indicates that u 2 (t) acquires a significant size as 
time goes on, see Fig. 3. 

Moreover, u(x,t) does converge to s (x) as t — >■ oo, contrary to the rather common 
belief that the most unstable Fourier more around u = (in our case, e iq2X ) will always win if 
it is excited. This example shows that the global evolution of the perturbation is not correctly 
described by the local analysis around the starting point, and that the final state is difficult 
to predict in general. However, under suitable assumptions on the parameters q , q 1 , these 
difficulties can be overcome and the existence of a trajectory connecting s qo to s can be 
effectively proved, as we now proceed to show. 

We begin by defining the function space in which we shall control the evolution of the 
system. According to Eq.(3.2), we let U qo s be the Hilbert space of functions u(x) of the form 

u(x) = uj^-^ , u n e C , (3.3) 

equipped with the scalar product 

1 f T 

i u , v )u = ™ (uv+u'v')dx = ^2(1 + (q - nSf)u n v n , 
1 Jo nez 

where T = 2it/8. It is easy to see that U qo s is invariant under the Ginzburg-Landau evolution 

u = u" + u — ul^l 2 , (3-4) 

and that u G U qo 6 if and only if e~ iqoX u(x) is periodic of period T and belongs to the Sobolev 
space H 1 ([0, T]). So, U qo ^ is isomorphic to the space 

Hj o ([0,T]) = {«GH 1 ([0,T])| 11 (T) = e i ^(0)} . 




Fig. 3: Time evolution of the three main Fourier modes u n (t) of u(x, t) in the case q Q = y/2/5, 8 = q Q /2. The 

periodic state s q (x) is slightly perturbed by ee qiX , where q l = q — 8, and the system undergoes a transition to 
the state s q ^{x). In the intermediate time lapse, the amplitude u 2 acquires a significant size, about 10% of the total 
amplitude. The dashed lines are the three projections of the solid curve onto the coordinate planes. 

Next, we note that the Ginzburg-Landau equation (3.4) is a "gradient system" in U qo s , in 
the sense that u = —8H/8u, where 

H(u) = ^J\\\u>?-\\u\> + \\uf) dx , T= f . 

In particular, the "energy" H (u) is a decreasing function of time when u is evolved with the 
GL equation. Our idea is to use this functional to show that any perturbation of s (x) in 
M q0: s which decreases the energy H will drive the system to the stationary solution s qi (x). In 
particular, any initial condition u(x,0) belonging to the unstable manifold of s qo (%) in U qo 6 
converges to s qi (x) ast — >■ oo. 

Forour argument to work, it is essential that s qi is the only stationary solution u ofEq.(3.4) 

in U qo 6 satisfying H{u) < H(s qo ). Thus, in addition to q > l/y/3 and 8 < «J 6<?q — 2, we 
have to assume that q < 8, for if q > 8, then the wave-number q 2 = q — 28 verifies \q 2 \ < q , 
so that the stationary solution s q2 (which belongs to U qo 6 ) satisfies H(s q2 ) < H(s ). Note, 
however, that the assumption q < 8 excludes some interesting cases, like the one represented 
in Fig. 3. 
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According to these remarks, the optimal parameter range within the scope of our method 

is 

£ = {(?„,*) I % < 8 < a/ 6? o " 2 and < % < 1} • 

This domain is represented in Fig. 4 in the q ,q 1 plane, where q 1 = q — 8. Although our 
results probably hold for all (q , 8) £ £ , they are somewhat easier to prove when 8 > 1, and for 
simplicity we shall restrict ourselves to this case. Thus, we define 



£ = {(q ,8) | 1 < 8 < J6ql - 2 and 1/^2 < q < l} C £ . 




Fig. 4: The regions of the q l , q plane corresponding to the domains £ (shaded) and £ C S (black). For all q , q l in 
the black region, we prove the existence of an Eckhaus transition from s to s , and we believe that the same is true 
for the whole shaded region. 



Proposition 3.1. If (q , 8) £ £ , then the only elements u in U qo 6 which are stationary solutions 
of the Ginzburg-Landau equation and whose energy H(u) satisfies H(u) < H(s qo ) are the 
functions 

u(x) = e^^l-qle^ , 

where q 1 = q — 8 and <p £ [0, 2tt). 

Using this result (which will be proved below), we now state and prove the main result of 
this section: 

Theorem 3.2. Let(q ,8) e £. Letu w e U qo S , and assume thatH(u w ) < H(s qo ). Then, as 
t — >■ oo the solution of the GL equation with initial data converges to 



u(x) = e lip y 1 — q\ e iqiX for some ip £ [0,27r) , 
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where q 1 = q — 8. 

Remark. Although it is linearly unstable, the fixed point s qo has an infinite-dimensional center- 
stable manifold V cs in U qo s , so there exist many perturbations of s qo in U qo s for which the 
conclusion of Theorem 3.2 fails. For example, (1 + e)s qo will converge to s qo , and not to 

s„ e l(p . On the other hand, there is a small neighborhood N of s„ in U„ c such that, for any 

9i ' 9o 9o,o ' J 



U 



^ £ N\V CS , the solution u t converges to s e l<p as t — >■ oo. Indeed, it is well-known from 
the theory of dynamical systems that such a trajectory leaves a (bigger) neighborhood of s qo 
essentially along the unstable manifold, and thus satisfies the energy condition H(u t ) < H(s qo ) 
for t sufficiently large. In particular, we can formulate the following 

Corollary 3.3. The conclusion of Theorem 3.2 holds ifu i0 \x) = s qo (x) + ee iqiX with e 
sufficiently small. 

Proof of Theorem 3.2. Denote by u t the (unique) solution of Eq.(3.4) in U qo s with initial data 

u,(°\ As a solution of a parabolic equation with smooth coefficients, u t is a smooth function of 
x,t for all t > 0. In particular, H(u t ) is a C 1 function of t, and an easy calculation shows that 

di H ( u t} = ~T J o \ u " + u t ~ \ u t\ 2u t\ 2 dx = -j! J o \u t \ 2 dx , 

where T = 2it/8. Since H(u) > —1/4 for all u £ M qoj s, we see that H(u t ) is decreasing and 
bounded from below, hence converges as t — >■ oo to some value . 

On the other hand, it is well-known that the trajectory {u t } t>0 is precompact in U qo s , see 
for example [Te]. In particular, its u;-limit set fi is non-empty, and any u £ fi is an equilibrium 
point of the system satisfying H(u) = < H(s ), see [Ha]. By Proposition 3.1, Q is thus 

contained in the circle of fixed points s qi e lip , <f £ [0, 2tt]. Now, a straightforward application 
of the stable manifold theorem shows that this circle has a small attractive neighborhood A in 
U qo s , and that every trajectory of the system entering A converges to a point of the circle as 
t — >■ oo [EG]. By definition of the a;-limit set, this is the case for the trajectory u t , and the proof 
of Theorem 3.2 is complete. 

Proof of Proposition 3.1. We first recall that the stationary solutions u(x) of the Ginzburg- 
Landau equation form a two-parameter family which can be indexed by the values of two "first 
integrals" J and E, see [DGJ, BR, Ga]. If u(x) is written in polar coordinates r(x)e lip ^ x \ then 
E and J can be expressed as 

J = r 2 (x)<p'(x) , E = ||ti'| 2 (a;) + ||ti| 2 (a;) — ||ti| 4 (a;) 

J 2 (3.5) 

1 If \ 2 I 1 i 1 2/ \ 1 4/ \ v ' 

= 2 r ( x ) + 2 + 2 r ( x ) ~ V ( x ) ■ 



Note that these formulas allow for an interpretation of E and J as the energy and the angular 
momentum of a point particle (in R 2 ) moving in the radially symmetric potential V(r) = 
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\r 2 — |r 4 . In this mechanical analogy, the role of the time is played by the space variable 
x £ R. 

As is well-known (see, e.g, [Ga]), the only stat ionary s olutions of Eq.(3.4) whose modulus 
|ii(a:)| is constant in x are the functions s (x) = y/l — q 2 e iqx for q G [—1, 1]. Since ^ > 1, it 
is easy to verify that s q £ U 6 if and only if q = q , q = q x , or q = 1. In this last case, s q = 

andO = H(s ) > — f (1 — 1o) 2 = H( s qo )- So, we see that Proposition 3.1 holds if we consider 
only stationary solutions with constant modulus. 

To conclude the proof, we need therefore only show that the space U qo s contains no 
stationary solution of Eq.(3.4) whose modulus is periodic in x and non-constant. If u(x) = 
r(x)e lip ^ is any such solution and if J, E are the corresponding values of the first integrals 
(3.5), then the (minimal) period of the modulus r(x) is given by 



T(E,J) = 2 [ 

J r 



^- J ) dr 

ri {E,J) (2E- J 2 r" 2 -r 2 + |r 4 ) 1/2 



(3.6) 



where r 1 and r 2 are the two smallest roots of the denominator, see [DGJ]. In the mechanical 
analogy above, this period corresponds to the time between two minima (or maxima) of the 
motion of the particle. Another important quantity is the total increase of the phase <p(x) during 
a period T. It is given by 



rr 2 (E,J) j 7 

viE, J) = 2 - — ^ , (3.7) 



and corresponds to the "advance of the perihelion" in the mechanical analogy. 

Remark. Strictly speaking, the formulas (3.6), (3.7) make sense only if J ^ 0. If J = 0, then 
u(x) is real (up to a global phase factor) and it is easy to verify that such a function cannot 
belong to U qo s if (q , 6) £ £ , except for the trivial case u = which we have already excluded. 
Thus, we shall always assume that J/0. 

It follows easily from the definitions that a stationary solution u(x) = r(x)e lip ^ of 
Eq.(3.4) parametrized by E, J belongs to U qo s if and only if 

jT(E,J) = ^ and jr](E, J) = mo d 2ir , (3.8) 

o o 

forsomej G N. In other words, the modulus r(x) has to be periodic of (not necessarily minimal) 
period 2it/8, and the increase of the phase <p(x) over a period must be equal to q (2ir/S) modulo 
27r. However, the next two lemmas will show that this is not possible if (q , 6) £ £ . 

Lemma 3.4. Let E,J be as above, and let < y 1 < y 2 < y 3 be the roots of the polynomial 

y 3 - 2y 2 + 4Ey - 2 J 2 = (y - Vl )(y - y 2 )(y - y 3 ) . (3.9) 

Then the period T(E, J) satisfies 

T > * » • ( 3 - 10 ) 

V 3 ^3 - 2 
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In particular, since y 1 + y 2 + y 3 = 2, one always has y 3 < 2, hence T > -n. 

Lemma 3.5. Let E, J as above, and assume that 8 = 2tt/T(E, J) > 1. Then the phase increase 
r](E,J) satisfies 

' < |'| < *VI + ^- (3 - n) 

In particular, one always has ir < \r)\ < ^/2-k. 

Using these two lemmas (which will be proved in the Appendix) we can now conclude 
the proof of Proposition 3.1. Let (q ,8) £ £, and assume that there exists a stationary solution 
u £ U 6 of the GL equation (3.4) whose modulus \u\ is periodic in x (and non-constant). If 
E, J are the corresponding values of the first integrals (3.5), then Eq.(3.8) is verified for some 
j £ N. Since 8 > 1 and (by Lemma 3.4) T > ir, the first equality in Eq.(3.8) can only be 
satisfied if j = 1. The second equality then implies that 77 = 2itq /8 + 2itm for some m £ Z, 
and by Lemma 3.5 we must have m = — 1 if 77 < and m = if 77 > 0. Combining these facts 
with the bounds (3.1 1), we obtain the inequalities 

f fi P fi ¥ 

In particular, we find 8 > ^Qql — 2, in contradiction with the definition of £. This means that 
Eq.(3.8) cannot be satisfied if (q , 8) £ £ , and the proof of Proposition 3.1 is complete. 

4. The Swift-Hohenberg Equation 

In this section, we address the question of phase slips in the Swift-Hohenberg equation 

d t U(x,t) = (3s 2 - (1 + d 2 x ) 2 )U(x,t) - U 3 (x,t) , (4.1) 

where U : R x R + — >■ R. For small e > 0, this equation has a family of spatially periodic 
stationary states of the form 

S = 2xcos(ux) + 0(x 2 ), (4.2) 

where 1 and 3e 2 = (lo 2 - l) 2 + 3 x 2 , see [CE1]. These states are known to be unstable 
if (u> 2 — l) 2 > e 2 and marginally stable otherwise. Numerical calculations show that, when 
such unstable states are perturbed, they typically relax by undergoing phase slips which change 
the number of zeros. A natural setup, which is considered in theory and experiments alike, 
consists in applying space periodic perturbations to these solutions. Then the phase-slips occur 
at periodically spaced points, changing in this way the density of zeros. In experiments this 
jump in density is in fact used to determine the limits of stability of the periodic solutions. 

The aim of this section is to present a rigorous description of these phenomena, using the 
results of the preceding sections for the GL equation. The relation between the two problems 
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is the so-called "approximation property," which says that the evolution of the SH equation is 
accurately described by the GL equation for a long (but finite) interval of time, see [CE2, vH, 
KSM, Sch]. A precise formulation of this property is as follows: First, fix T > 0, R > 0, and 
let t) be a solution of the GL equation 

d T A(Z,T) = dlA(^r) + A(£,r) - A(^t)\A(^t)\ 2 , (4.3) 

satisfying || 0)|| c4 < R. Next, for any e > 0, define 

v a,A x ^) = 2eRe (A(^ex,Ze 2 t)e ix ^j > ( 4 - 4 ) 
and let U A (x, t) be the unique solution of the SH equation in C 4 (R) with initial condition 

U A ,A^) = V AtS (x,0). (4.5) 

Under these assumptions, there exist an e > and a C < oo such that, for all e < e and all 
< G [0,T /e 2 ], one has 

< c= 4 *. ( 4 - 6 ) 

see [KSM]. In particular, the approximation V A£ (-,<) is within C(e 2 ) of the true solution 
U Ate (-,t) when^G [0,T /e 2 ]. 

Remark. In [KSM] , the authors consider the more general case where U Ae (x,Q) = V A e (x,Q) + 
0(e 2 ), and they show that \\U A ^(-,t) - V A ^(-,t)\\ < Ce 2 for all t e [0,T /e 2 ]. However, 
in the particular case where the initial data for U A ^ S and V A coincide, their proof yields the 
stronger result (4.6). 

We first give a local description of a typical phase slip for the SH equation, using the results 
of Section 2. Let A(£,r) be a solution of the GL equation (4.3) satisfying the assumptions of 
Theorem 2.7 for u(x,t). Without loss of generality, we may assume that A(0, d) = 0, i.e. the 
phase slip occurs exactly at £* = 0, r* = d. Then, according to Eq.(2.3), there exist a > 0, 
b £ C, Re b ^ 0, and d > such that 

A(£, r) = \ai 2 + ibi + a (r-d) + (|£| 3 + \r - d\ 3 ' 2 ^ , (4.7) 

for all \i\ < d 1/2 and all r e [0, 2c?]. For e > 0, we define V A s {x,t) by Eq.(4.4) and we let 
U A e (x,t) be the solution of the SH equation satisfying the initial condition (4.5). If d < 3T /2 
and e is sufficiently small, then Eq.(4.6) implies that 

V A , e (x,t) = V Ate (x,t) + O(eH) , (4.8) 

for all \x\ < 2(d/3) 1/2 e~ 1 and all* G [0,2<Z/(3e 2 )]. 

To show that U(x,t) undergoes a phase slip near x = 0,t = d/(Se 2 ), we assume from 
now on that e < d/K for some (large) constant K > 0, and we restrict Eq.(4.8) to the smaller 
region 

R d,e = {OM) I l«l < -K" > \t~ d/{Se 2 )\ < K/e) . 
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Setting b' = b[ + ib' 2 = Vsb/2, we find 

U As (x,t) = 2e 2 Re ((ib'x + a(Set - d/e))e ix ) + G(e 2 d) 

= 2e 2 ((a(3et — d/e) — & 2 a:)cos x — b\ x sin x) + 0(e 2 d) , 

uniformly in R d e . 

Since we wish to study how the number of nodes of the solution changes with time, we 
look for points where U A e (x,t) = 0. The leading term in Eq.(4.9) vanishes if and only if 
f b (x) = a(Set — d/e), where 

f b (x) = b\x tan x + b' 2 x . (4-10) 

The graph of the function f b (x) — s is represented in Fig. 5 for b\ < 0. We see that if —s ^> 1, 
the roots of the equation f b (x) = s are located near the points x = (n + l/2)7r, n £ Z. As 
s goes through zero, the two central roots coalesce and disappear in a neighborhood of the 
origin. Conversely, if b[ > 0, a pair of roots is created near x = 0. These phenomena are 
precisely what we call phase slips for the SH equation. They will occur if we choose the constant 
K = K(a, b) > so that the quantity s = a(3et — d/e) varies over a sufficiently large range 
as t varies in the region R d s . Of course, when considering the full Eq.(4.9), the picture will 
be slightly distorted by the error terms, but qualitatively it remains the same if d is sufficiently 
small. These arguments provide a sketch of the proof of the following 



fb( x ) ~ s 




:> x 



Fig. 5: The amplitude of the SH equation, divided by cos a;. We show the graph of the function f b (x) — s = 
b\ x tan x + bLx — s for s < 0. 
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Theorem 4.1. Assume that A(£,t) satisfies the assumptions of Theorem 2.7 for u(x,t). For 
e > 0, let U A (x,t) be the solution of the SH equation with initial condition 

2eRe (a(3 1/2 sx/2, 0)e ix ^j . 

Then there exist a d 1 < d* and a K = K(a,b) > such that, for all d £ [0,^] and all 
e < d/K, the function U A e (x,t) undergoes a phase slip near x = 0, t = d/(3e 2 ). This phase 
slip corresponds to the creation (Re b > 0) or the annihilation (Re b < 0) of a pair of roots of 

We next address the question of the global evolution of an Eckhaus instability for the 
Swift- Hohenberg equation. Assume that q ,q 1 are two wave-numbers such that (q ,q — ?i) 
belongs to the region £ defined in Section 3, and let Sj(£) = (1 — q^) 1 ^ 7 'e iqj ^ for j = 0, 1. 
If e > is sufficiently small, the SH equation (4.1) has two families of periodic stationary 
solutions Sj (x) which are of the form 

S 3 Jx) = 2eRe (s 3 {4exY^) + 0(e 2 ) . (4.11) 

We recall that S (x) is always linearly unstable, whereas S 1 ^(x) is stable or unstable ac- 
cording to whether q\ < 1/3 or q\ > 1/3. Combining Theorem 3.2 with the approximation 
property (4.6), we shall show the following 

Theorem 4.2. Assume that (q , q — q x ) e £ . For any p > 0, there exist e > 0, T > such 
that, for all e £ [0, e ], the Swift-Hohenberg equation (4.1) has a solution U(x,t) satisfying 

\\U(;0) ~ S 0|0 (-)|| L oo < ep , \\U(.,Tjs 2 ) - S.J^^ < ep , (4.12) 

for some <p G [0,27r). 

Remark. The construction of U(x,t) will show that U(x, 0) and S (x) have the same density 
of zeros (if p is sufficiently small), and similarly for U(x, T /e 2 ) and S 1 (x). Since q ^ q 1 , 
this means that U(x,t) undergoes phase slips as t varies from to T /e 2 . 

Remark. Since the solutions depend continuously on time, there is in fact a neighborhood of 
initial data near U(-,0) for which the conclusions of Theorem 4.2 hold. In this sense, the 
Eckhaus instability leads generically to phase slips. 

Proof. Fix p > 0. According to Theorem 3.2, there exists a solution A(£, r) of the GL equation 
(4.3) such that — oo) = s (£) andA(£,oo) = e l ' p s 1 (£), for some <p £ [0,2-7r). Therefore, 
there is a T > such that, after choosing the origin of r conveniently, one has 

- -o(0IIl<- < P/ 3 , \\M;3T ) - e^ Sl {.)\\ Lx < p/S . (4.13) 

Furthermore, by the approximation property, there is an £ > such that, if V(x,t) = V A (x,t) 
is defined by Eq.(4.4) and if U(x,t) = U A (x, t) is the solution of the SH equation with initial 
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condition (4.5), then Eq.(4.6) holds for all e < e . In view of Eqs.(4.6), (4.11), (4.13), we thus 
have 

l|tfM)-S 0|0 (.)|| £ «» = |in-,0)-5 0|0 (.)|| £ «- 

< 2e||A(-,0) -s (-)\\ Loo + 0{e 2 ) < 2ep/3 + C*e 2 , 

and similarly 

\\U(;Tje 2 ) - S ltV (.)\\ Loo < \\U(;Tje 2 ) - V{-,TJ s 2 )\\ Loo + ||V(-,r o / £ 2 ) - S ltV (.)\\ La 

< Ce 2 + 2e\\A(.,3T )-e i *s 1 (.)\\ Loo + (D(e 2 ) 

< C*e 2 + 2ep/S . 

Therefore, if C*e < p/3, the solution U(x,t) verifies Eq.(4.12). The proof of Theorem 4.2 is 
complete. 

Remark. Unfortunately, it is not possible to conclude from the preceding arguments that the 
solution U(x,t) actually converges to S 1 (x) as t — >■ oo. Indeed, the approximation (4.6) is 
only valid over a finite interval of time, whereas the trajectory r) takes an infinite time to 
reach the fixed point e l(p s Moreover, even when q\ < 1/3, we do not know if the stationary 
solution S 1 (x) of the SH equation is non-linearly stable. 



Appendix 



Proof of Lemma 3.4. We first rewrite the expression for the period T in a more convenient 
form. Setting y = r 2 in Eq.(3.6), we find 

T _ f V2 dy _ ^ f y2 dy 

~ J yi y/2Ey-J* -y>+y*/2 ~ 2 J yi y/(y - yi )(y 2 y)(y 3 - y) ' 

Using the change of variables 



V = a (V0 = y x cos 2 (V>) + y 2 sin 2 (V>) , dy = 2 yj (y - y 1 )(y 2 - y) dip , 
we obtain 

T = 2^' r/2 ** ^ 



Jo 



^/y 3 - a(V>) M(y/y 3 - y x ^y z ~ Vz) ' 



where A4(a, b) is the arithmetic-geometric mean of a and b, see [BB]. By the definition of this 
mean, we have M(a, b) < (a + b)/2, so that 

T > - . A.l 

V^s - yi + V^s - ^2 



Eckhaus Instability 



20 



It remains to show that Eq.(A.l) implies Eq.(3.10). The idea is to fix y 3 and to consider 
the right-hand side of Eq.(A.l) as a function of y 1 alone, using the fact that y 2 = 2 — y 1 — y 3 . 
Thus, we define 



F (yi) = Vy^ -vi + Vy^ - y-i = -vi + \J 2 v% - 2 + vi > 

and an easy calculation shows that F\y 1 ) > 0. Now, since y 1 < y 2 and y 1 + y 2 + y 3 = 2, we 
always have y 1 < 1 — y 3 /2, hence F(y 1 ) < F(l — y 3 /2) = -\J^y 3 — 4. Inserting this bound 
in Eq.(A.l), we find 

271-^2 2?r 
as asserted. The proof of Lemma 3.4 is complete. 

Proof of Lemma 3.5. Since the GL equation (3.4) is invariant under the complex conjugation 
u — > u, there is no loss of generality in assuming that J > 0, hence 77 > 0. Proceeding as in the 
proof of Lemma 3.4, we rewrite the expression (3.7) of 77 as 

77 = 2V2 / — — . Y , (A.2) 
Jo a(V>) ^3 - «W 

where again a(ip) = y 1 cos 2 (ip) + y 2 sin 2 (ip). 

To prove the lower bound on 77, we first note that y-yy 2 y 3 = 2 J 2 , in view of Eq.(3.9). On 
the other hand, an elementary calculation gives 



2 r /2 dij> 



Using Eq.(A.2), we obtain 

r / 2 # 



(A.3) 



»7 > 2 Vy^y^h 

Jo 



7T 



To prove the upper bound on 77, we apply the Schwarz inequality to the right-hand side of 
Eq.(A.2) and obtain the bound 

1/2 '- -*/. ^ V /2 



77 7I-V2J ( ^ ^ j ( TT / 



\t Jo o(^) y \ir Jo a(i>)(y 3 - a (V0) J 

Using Eq.(A.3) and the analogous result 

2 r /2 # _ j_ / 1 1 \ 

Wo o(^)(y 3 ~ a (V0) j/s w^i^ V^s -vWva 
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we thus find 

"^(W (»-X-» tT- (a - 4) 

It remains to show that Eq.(A.4) implies the upper bound in Eq.(3.1 1). First, we note that 
Eq.(3.10) implies ■ S /Sy 3 — 2 > 2-k/T = 8, hence y 3 > (2 + 8 2 )/S. Since 8 > 1 by assumption, 
we see that y 3 > 1 too. Next, proceeding as above, we fix y 3 £ (1,2) and define 



(y 3 - yi)(y 3 - y 2 ) (y 3 - 2/i)( 2 2/s - 2 + Vi) ' 

By direct calculation, we find 

n i ( ^ 2 (y 2 -yi)y 3 (y 3 - 1) . n 

G(Sl)= ( S 3- Sl F( S 3- S2 ) 2 

so that G(y 1 ) < G(l — y 3 /2) = (2 — y 3 ) 2 /(Sy 3 — 2) 2 . Since this expression is a decreasing 
function of y 3 , we can replace y 3 by the lower bound (2 + 8 2 )/S. We thus obtain 

G ^ (2/3 — 2/i)(2/s — 2/2) ( 3 ^ 2 ) ' ^ A ' 5 ' ) 
Combining Eqs.(A.4), (A.5), we find 



4 - 8 2 2 4 

?7<7r V 1 + ^ = 7r V3 + ^' 

as asserted. The proof of Lemma 3.5 and hence of the global picture for the GL equation are 
complete. 
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